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Abstract 

A study of differential event shape distributions using e + e~ data at centre-of-mass 
energies of ^fs = 35 to 183 GeV is presented. We investigated non-perturbative power 
corrections for the thrust, C-parameter, total and wide jet broadening observables. 
We observe a good description of the distributions by the combined resummed QCD 
calculations plus power corrections from the dispersive approach. The single non- 
perturbative parameter cto is measured to be 

a (2 GeV) = 0.502 ± 0.013(stat.) t[J;^(exp. syst.) ±g;g£f(theo. syst.) 

and is found to be universal for the observables studied within the given system- 
atic uncertainties. Using revised calculations of the power corrections for the jet 
broadening variables, improved consistency of the individual fit results is obtained. 
Agreement is also found with results extracted from the mean values of event shape 
distributions. 
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1 Introduction 

The study of final states in electron-positron annihilation led to the development of numer- 
ous observables as powerful tools for the analysis of hadronic events. Improving theoretical 
predictions for these observables in the last years as well as the large amount of data col- 
lected by various experiments e.g. at the PETRA, PEP, LEP and SLC accelerators from 
^Js = 12 to 189 GeV provide precise quantitative tests of Quantum Chromodynamics 
(QCD). The key feature of perturbative QCD, the running of the strong coupling constant 
a s , has meanwhile become an indubitable experimental fact from e + e~-data (see [1] and 
references therein), in particular due to investigation of the energy dependence of event 
shape data. 

Precision tests of perturbative QCD from hadronic event shapes require a solid un- 
derstanding of non-perturbative effects. In the context of hadronisation, non-perturbative 
effects are usually estimated by phenomenological hadronisation models available from 
Monte Carlo event generators. Alternatively, analytical approaches were pursued in the 
past years in order to deduce as much information as possible about hadronisation from the 
perturbative theory, e.g. the renormalon inspired ansatz [2] or the dispersive approach [3]. 

Phenomenological considerations lead to the expectation that hadronisation contribu- 
tions to event shape observables evolve like reciprocal powers of the hard interaction scale 
^/s. The investigation of power corrections therefore substantially benefits from e + e~data 
at the lower energy region as marked by PETRA/PEP experiments at ^/s = 12 to 47 
GeV. Our re-analysis of e + e~data collected by the JADE detector [4,5], one of the former 
PETRA experiments, allowed studies of the energy evolution of QCD at combined lower 
and higher energies in terms of observables frequently used in the analysis of events at the 
highest centre-of-mass energies available at LEP. Further efforts were undertaken to make 
event shapes available at hadronic centre-of-mass energies below the Z pole, e.g. the anal- 
ysis of radiative events recorded at LEP-1 as performed by several LEP experiments [6]. 

In the present paper we test power corrections to the differential distributions of event 
shape observables measured from re-analysed JADE data at 35 to 44 GeV and from data 
of the LEP/SLC experiments at 91 to 183 GeV. We focussed on two-loop predictions [7,8] 
basing on the dispersive approach. For the jet broadening measures, progress was made 
in the understanding of non-perturbative corrections to the distributions [8], in particular 
the interdependence between perturbative and non-perturbative contributions. We also 
test a basic concept of the power corrections, the universality of the zeroth moment of the 
effective coupling a below a given infrared energy scale ///. 

Section 2 starts with an overview of the observables and the data used in the analysis. 
After some annotations to the parametrisation of the power corrections in Section 3, 
we present fit results basing on combined power corrections plus perturbative QCD in 
Section 4 and draw conclusions from our results in Section 5. 



2 Event shape data 

For our studies, we considered the event shape distributions of thrust, the C-parameter, 
the total and wide jet broadening, Bt and By/- For convenience, we list the definitions of 



these observables. 

Thrust T: 

The thrust value is given by the expression [9] 

rp [Ei\Pi-n\ 

1 = max 
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The thrust axis Ht is the vector n which maximises the expression in parentheses. 
It is used to divide an event into two hemispheres H\ and H 2 by a plane through 
the origin and perpendicular to the thrust axis. 

Jet Broadening B: 

The jet broadening measures are calculated by [10]: 

B = (EieH k \PiXnT\\ 



for each of the two hemispheres, H^, defined above. The total jet broadening is given 
B T — B 1 + B 2 . The wide jet broadening is defined by B w = max(5 1 , B 2 ). 

C-parameter: 

The C-parameter is defined as [11] 

C = 3(AiA2 + A2A3 + A3A1) 
where A 7 , 7=1,2, 3, are the eigenvalues of the momentum tensor 
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Experimental data for the thrust below the Z pole are provided e.g. by the experiments 
of the PETRA (12 to 47 GeV), PEP (29 GeV) and TRISTAN (55 to 58 GeV) colliders. 
For the thrust, we considered measured distributions from the PETRA experiments JADE 
and TASSO (total of about 50000 events), from the PEP experiments MARKII and HRS 
(total of about 28000 events) and from the AMY detector (about 1200 events) at TRISTAN 
(see Ref. [12], and [4]). For the other observables, we used experimental distributions at 
PETRA energies from our re-analysis of JADE data [4, 5] providing about 22000 events at 
35 GeV and 6200 events at 44 GeV. The distributions used at the Z resonance stem from a 
measurement using several hundred thousands events for each of the four LEP experiments 
OPAL, ALEPH, DELPHI and L3, and of about 40000 events for the SLD experiment at 
SLAC [13]. We also considered measurements performed by the LEP collaborations above 
the Z pole [14] providing a few thousend events per experiment up to a center-of-mass 
energy of 183 GeV. 

All distributions used in this study were corrected for the limited resolution and ac- 
ceptance of the respective detectors and of the event selection criteria (for details see the 
referenced publications). 



3 Power corrections to differential distributions 

Non-perturbative effects to event shape observables arise due to the contributions of very 
low energetic gluons which cannot be described perturbatively because of the presence of 
an unphysical Landau pole in the perturbative expression of the running coupling a s (Mzo). 
Within the dispersive approach [3, 7] a non-perturbative parameter 

1 rvi 
Oio(fJ-i) = — d/c a s (k) 
Hi Jo 

is introduced to parametrize the analytic but unknown behaviour of a s (y/s) below a certain 
infrared matching scale fij. The technique presented in the cited papers also removes the 
divergent soft gluon contributions in the perturbative predictions for the observables and 
produces as a consequence 1/y/s corrections. 

The power corrections which we used in this analysis have been calculated in Ref. [7,8] 
up to two loops for the differential distributions of the event shapes considered here. The 
effect of hadronisation on the distribution obtained from experimental data is described 
by a shift of the perturbative prediction away from the kinematical 2-jet region, 

da cxp dcr pert 

dT y ' dT y ' K ' 

where T = 1 — T, C, Bt and B\y, respectively. The factor V containing the only non- 
perturbative parameter a is the same for all observables, 

V = ^M^(aM-a s (^s)-Pof(ln^ + ?- + l)) (2) 

7T 2 y/S \ 27T \ 1*1 PO J J 

where Cp = 4/3. (3q = (HCa — 2iV/)/3 stems from the QCD renormalisation group 
equation for Nc = Ca colours and Nf active quark flavours. K = (67/18 — tt 2 )Ca — 5/9Nf 
originates from the choice of the MS renormalisation scheme. The Milan factor J\4 accounts 
for two-loop effects. Its numerical value of 1.795 (for three active quark flavours) is given 
by the authors of [7] with a theoretical uncertainty of about 20% owing to missing higher 
order corrections. Djr is a function depending on the observable, 

{2 ;T= 1-T 

3tt -T = C ( 3 ) 

\\n^ + B T {T,a s {T^)) ;F = B T ,B W . 

Thus a simple shift is expected for thrust and C-parameter, whereas for the jet broaden- 
ing variables an additional squeeze of the distribution is predicted. Besides the logarithmic 
dependence on the value T of the observable, the term B? exhibits an additional (weaker) 
dependence on the value of the observable and on the running coupling at the scale T\fs. 
The more complex behaviour for the jet broadening recently calculated in [8] is related 
to the interdependence of non-perturbative and perturbative effects which cannot be ne- 
glected for these observables. As a consequence, the term VDjr leads not only to a larger 
shift but also to a stronger squeeze at decreasing centre-of-mass energies. The necessity 
of an additional non-perturbative squeeze of the jet broadening distributions was already 
pointed out in [15]. 



4 Tests of power corrections 

A few experimental tests of power corrections to differential distributions have been done 
up to now [8,15-18]. Some studies for the new observables which were performed using 
LEP data only suffer from the low data statistics above the Z pole and from the lack of 
data at energies below the Z pole, where QCD scaling violations and energy dependent 
non-perturbative effects are larger. The present analysis of combined LEP and PETRA 
data allows a more stringent test of the 1/y/s dependence of the power corrections. 



4.1 Fits to the data 

The standard analysis uses the entire data set described in Section 2. For the perturba- 
tive predictions, we employed combined second order [19] plus resummed [20-23] QCD 
calculations (0(a 2 s ) + NLLA). Predictions of this type take into account the dominating 
leading and next-to-leading logarithms to all orders in a s in the perturbation series for the 
observable and therefore exhibit a weaker dependence on the QCD renormalisation scale \x 
= x^^. Since the matching procedure of the exact matrix elements and the NLLA calcu- 
lations reveals some ambiguities, we apply the four different matching schemes proposed 
in [22], the R-, ln(R), the modified R and the modified In (R)- matching. 

For each observable, we performed simultaneous x 2_ fits of both perturbative predic- 
tions plus power corrections which are parametrised by a s (M z o) and ao(fJ>i), respectively. 
a s (M z o) was evaluated to the appropiate energy scale ji = x^\/s using the two-loop for- 
mula for the running coupling [24]. The errors used in the fit are the quadratic sum of 
statistical and experimental systematic uncertainties of the event shape distributions de- 
noted in the references. The errors of the fit are derived from the 95% CL two-dimensional 
contour. The fit ranges are defined individually for each center-of-mass energy such that 
they, in general, include the maximum of the differential cross section located in the 2-jet- 
region of the distributions. In doing so we ensure to exploit the range of the distribution 
at each energy point as far as possible in order to constrain the theory and to consider as 
much data statistics as possible from the distributions, which is in particular important 
in regard of the scarce LEP-2 data. 

Our standard result for an observable is defined by the unweighted average of the 
results obtained from the different (!)(«;?) +NLLA-combination schemes. We fixed the 
renormalisation scale x^ — 1 and the infrared matching scale \xj — 2 GeV. The fit error 
assigned is the average over the individual fit errors for the matching schemes. The numeric 
results are listed in Table 1 for a s and in Table 2 for a . The fit curves for the modified 
ln(R)-matching and the corresponding experimental data for the thrust 1 — T, the C- 
parameter, the total and the wide jet broadening, B? and Bw, are shown in Figures 1, 
2 and 3. Generally we observe a good agreement of the predictions with the data within 
the fit ranges, independent of the matching scheme used in the fit. The agreement is 
moderate for some distributions at y/s < 35 GeV. The xVd.o.f . of the fits (see Table 3) 
are between 0.93 (for Bt, mod. ln(R)-matching) and 1.34 (for B\y, R-matching). The fit 
quality deteriorates for the jet broadenings if the upper limit of the fit range is extended 
too far into the 3-jet region. The excess of the theory over the data observed here has also 
been noticed in other studies [25-27] where event shape distributions were corrected for 



hadronisation effects using Monte Carlo models. Thus these discrepancies might be due 
to the perturbative predictions. 

The results for a s (M z o) obtained from the fits proved to be systematically lower than 
the respective results at y/s = M z o [5,25-28] which use the same perturbative predictions 
but apply Monte Carlo corrections instead of power corrections. They are consistent with 
each other within the statistical errors in the case of thrust, C-parameter and the total 
jet broadening Bj-, whereas the results obtained for By/ is about 20% below the world 
average value for a s (M z o) [1]. This deviation is to a minor extend also present in the 
studies cited. Nevertheless, using the improved power corrections the values for a s (M z o) 
for Bt and By/ are now much more compatible with the values obtained from the other 
observables than is the case in our previous study [15]. 

The numeric results for ao(2 GeV) vary between 0.437 for the C-parameter and 0.685 
for B\y. The high value for By/ is supposed to be partially related to the systematically 
lower values of a s , since the two fit parameters are anticorrelated (see Table 3 for the 
correlations coefficients) . 

From the experimental point of view, there is a lucid explanation for the deviations 
of the a s results based on the power corrections from those based on Monte Carlo cor- 
rections [15]. We realised in our studies that the latter induce a stronger squeeze to all 
distributions than the power corrections which for instance predict merely a shift for the 
thrust and the C-paramter without any presence of a squeeze. Although the situation 
was a little remedied for the jet broadenings B due to the factor In 1/B appearing in 
Equation (3), the extend of the squeezing remains below the expectation of hadronisation 
models. As a consequence, the two-parameter fit favours systematically smaller values for 
a s in order to make the predicted shape more peaked in the 2 jet region of the distribution 
(thus causing a squeeze) and hence chooses high values for a in order to compensate the 
resulting shift of the distribution towards the 2-jet-region, as is in particular the case for 
the jet broadenings. 

4.2 Systematic uncertainties 

In order to scrutinize the universality of the non-perturbative parameter a we also studied 
the impact of systematic uncertainties to our final results. Theoretical uncertainties from 
perturbation theory were assessed by considering the different matching schemes and by 
varying x^ from 0.5 to 2.0. For the further systematic checks, we took the modified ln(R)- 
matching scheme as a reference since it turned out that most of the relative changes of a s 
and a are n °t significantly affected by the choice of the matching scheme. Uncertainties 
due to the parametrisation of the power corrections are associated with the arbitrariness 
of the choice of the value of \xi which has been varied by ± 1 GeV. Furthermore the Milan 
factor Ai was varied within the quoted uncertainty of 20%. No error contribution from 
the infrared matching scale is assigned to a (/ i /)- 

The signed values in Tables 1 and 2 indicate the direction in which a s (M z o) and «o 
changed with respect to the standard analysis. It turns out that the non-perturbative 
uncertainties for a s (M z o) are negligible for each observable. This is in accordance with 
the expectation since a s is mainly constrained by the perturbative prediction and not 
by the non-perturbative shift oc fiiA4(ao — a s + 0{a 2 s )). On the other hand, the strong 



dependence of «o on M. is due to the obvious anticorrelation. In regard to the remarks 
made in Section 4.1 it is noticable that the value for a from Byy is affected with the 
largest non-perturbative error. 

We also examined the dependence of the results from the input data taken for the fits 
which is for brevity denoted as "experimental uncertainty". Since the statistical weight 
of the event shape distributions at ^/s = 91 GeV is dominant in the fits, the analysis 
was repeated using the data from each LEP/SLC experiment seperately as representative 
for y/~s = 91 GeV. Each deviation from the standard result was added in quadrature. In 
addition, we performed a measurement omiting the 91 GeV data. The larger of the two 
latter deviations were incorporated in the experimental error. Furthermore, an analysis 
was performed using seperately either event shape data below the Z pole or data at ^/i > 
Mzo. This checks also for higher order non-perturbative contributions oc (1/^/i) 2 to the 
power corrections. 

A further source of experimental uncertainty comes from the choice of the fit range. We 
quantified it by varying the lower and the upper edge of the fit range of all distributions 
likewise in both directions by an amount approximately corresponding to the respective 
bin of the JADE distribution at 35 GeV. The fit range of a distribution was expanded 
into the 2-jet region as far as the x 2 of the extreme bin did not contribute substantially 
to the total \ 2 of the distribution. For each edge, we took the larger deviation from the 
standard result as error. 

In general, the experimental errors are more moderate than the theoretical uncertain- 
ties. We realise major variations of a s (Mzo) from Bt and By/ when considering only 
PETRA data in the fits. The larger dependence of the results on the lower edge of the 
fit range of the thrust distribution may be partially ascribed to the contributions of the 
lower energy data down to 14 GeV at which mass effects become important. 

All errors are treated as asymmetric uncertainties on a s (Mzo) and ao{fii), respectively. 
The total errors are quoted as the quadratic sum of the individual uncertainties, thus 
disregarding possible correlations between them. They are also given in Table 1 and 2 for 
a s and a , respectively. 

4.3 Combination of individual results 

The individual results for the four event shape observables were combined to a single value 
for a s (M z o) and ao(2 GeV), respectively, following the procedure described in [4]. This 
procedure accounts for correlations of the systematic errors. A weighted average of the 
individual results was calculated with the square of the reciprocal total errors used as the 
weight. For each of the systematic checks, the mean values for a s (M z o) and «o(2 GeV) 
from all considered observables were determined. Any deviation from the weighted average 
of the main result was taken as a systematic uncertainty. 

With this procedure, we obtained as final result for ao(2 GeV) 

a (2 GeV) = 0.502 ± 0.013(stat.) la&xP- s y st -o:o53( theo - s y st 

where the error is dominated by non-perturbative uncertainties. Figure 4 illustrates that 
the individual results are compatible with the weighted mean within the total errors. We 



consider this result as a confirmation of the predicted universality of the non-perturbative 
parameter a . 

Employing the same procedure for a s (M z o), we obtain 

a s (M z o) = 0.1068 ±0.0011(stat.) l™j(exp. syst.)^;S(theo. syst.) . 

This small value compared to the world average [1] is due to the a s from Bw which is 
substantially below the a s from T, C-parameter Bt (see Figure 4). 

If B w is omited, the final value is a s (M z o) = 0.1141 ±0.0012(stat.) ^[^(exp. s y st 
-ooo4i(t neo - syst.) which is in good agreement with the world average value [1]. In this 
case, the result for a does insignificantly change to a (2 GeV) = 0.478 ±0.013(stat.) 
t^oVp. syst.) lg:8S(theo. syst.). 

5 Summary and conclusions 

The analytic treatment of non-perturbative effects to hadronic event shapes based on power 
corrections were examined. We tested predictions [7, 8] for the differential distributions of 
thrust, C-parameter, total and wide jet broadening, Bt and Bw, respectively. For this 
test, a large amount of event shape data collected by several experiments over a huge 
range of e + e~ annihilation energies from ^Js = 14 to 183 GeV was considered. For the 
C-parameter, Bt and Bw, differential distributions below the Z pole are provided by our 
re-analysis of e + e~data of the JADE experiment [4,5] at ^/i = 35 and 44 GeV. 

A simultaneous two-parameter fit of the combined power corrections plus perturbative 
QCD-predictions ((9(c^)+NLLA) was performed for both the strong coupling a s (M z o) 
and the non-perturbative parameter ao (/•*/)• The good quality of the fits with x 2 /d.o.f. 
around 1.0 for each observable supports the predicted 1/^/s-evolution of the power cor- 
rections. The results for a s (M z o) and «o(2 GeV) are more consistent with each other 
due to the improved predictions of power corrections to the jet broadening variables [8]. 
Nevertheless we still observe a large deviation of the results obtained from Bw from those 
extracted from the other observables. We suspect it to be partially due to the resummed 
predictions for this observable. The individual results for a s from all observables (Table 1) 
are systematically smaller than the respective results in [5,25-28], the latter employing 
hadronisation models for the correction of non-perturbative effects. This observation is 
presumably related to the different extends of the non-perturbative squeeze of the distri- 
butions predicted by the both ansatzes. We obtain as combined result 

a s (M z o) = 0.107 ±8:g8i . 

It should be noted that a s from the wide jet broadening Bw is rather small compared with 
the a s from the other observables. The combined value for a s (M z o) derived when omiting 
the results for Bw is a s (M z o) = 0.114 loioosO'Ot.), in agreement with direct measurements 
at the Z resonance [25-28] and with the world average value of a^ a '(M z o) = 0.119±0.004. 
The combination of the individual results for a from all observables considered yields 
as final result for the coupling moment 

a (2 GeV) = 0.50 ±°;°° 



where the error is dominated by theoretical uncertainties defined in Section 4.2. Since this 
value is representative for the individual results within the errors defined, we consider this 
as a confirmation of the universality of ao as predicted within the dispersive approach [7] 
of power corrections. Nevertheless, the results yielded from Bt and By/ indicate that 
higher orders may contribute significantly to the non-perturbative corrections for the jet 
broadenings. 

We point out that the average value for ao(2 GeV) derived from the differential event 
shape distributions is in good agreement with the result of ao(2 GeV) = 0.47 to 04 provided 
by our previous study [5] of power corrections to the mean values of event shapes. 
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Tables 





1-T 


c 


Bx 


B w 


a s (M z o) 


0.1156 


0.1137 


0.1125 


0.0973 


Exp. stat. 


± 0.0011 


± 0.0011 


± 0.0014 


± 0.0009 


pQCD 


+0.0055 
-0.0043 


+0.0051 
-0.0036 


+0.0063 
-0.0054 


+0.0026 
-0.0015 


Pow. corr. 


+0.0007 
-0.0006 


+0.0000 
-0.0001 


+0.0002 
-0.0003 


+0.0000 
-0.0000 


Exp. syst. 


+0.0034 
-0.0021 


+0.0026 
-0.0022 


+0.0050 
-0.0023 


+0.0015 
-0.0043 


Total 


+0.0066 
-0.0049 


+0.0058 
-0.0044 


+0.0082 
-0.0061 


+0.0031 
-0.0046 


ln(R) 


-0.0002 


-0.0001 


-0.0015 


< 0.0001 


ln(R) mod. 


+0.0008 


+0.0011 


+0.0029 


+0.0013 


R 


-0.0010 


-0.0005 


-0.0032 


-0.0010 


R mod. 


+0.0004 


-0.0005 


+0.0018 


-0.0003 


x^ = 0.5 


-0.0041 


-0.0035 


-0.0041 


-0.0010 


x^ = 2.0 


+0.0055 


+0.0050 


+0.0053 


+0.0022 


M - 20% 


+0.0003 


< 0.0001 


+0.0002 


< 0.0001 


M + 20% 


-0.0003 


< 0.0001 


-0.0001 


< 0.0001 


91 GcV = OPAL 


-0.0006 


+0.0004 


-0.0011 


-0.0023 


91 GeV = ALEPH 


-0.0001 


-0.0004 


— 


— 


91 GeV = DELPHI 


-0.0001 


+0.0005 


+0.0009 


+0.0011 


91 GeV = L3 


-0.0010 


-0.0011 


— 


— 


91 GeV = SLD 


-0.0014 


-0.0017 


+0.0004 


-0.0008 


y/s ^ 91 GeV 


-0.0003 


-0.0004 


+0.0020 


-0.0012 


y/s > 91 GeV 


+0.0019 


< 0.0001 


+0.0002 


+0.0001 


y/s < 91 GeV 


+0.0007 


+0.0024 


+0.0046 


-0.0035 


fit range 


+0.0027 
-0.0010 


+0.0004 
-0.0009 


+0.0006 
-0.0021 


+0.0010 
-0.0005 



Table 1: Values of a s (Mzo) derived using combined power corrections plus resummed 
QCD predictions to the event shape distributions of trust T, C-parameter, total and 
wide jet broadening Bt and By/- The standard result for each variable is defined as the 
unweighted average of the results from the different matching schemes (ln(R), ln(R) mod., 
R, R mod.) using x^ = 1 and \xj = 2 GeV and considering the entire data set. In addition, 
the statistical and systematic uncertainties are given. Signed values indicate the direction 
in which a s (M z o) changed with respect to the standard analysis. Theoretical uncertainties 
are given by the matching scheme ambiguity, the choice of rr M and fij and the uncertainty 
of the Milan factor M.. Experimental uncertainties are taken into account by varying 
the input data of the fit (see text). All errors are treated as asymmetric uncertainties on 
« S (M Z »). 



11 





1-T 


c 


Bx 


B w 


a 


0.469 


0.437 


0.562 


0.685 


Exp. stat. 


± 0.011 


± 0.010 


± 0.017 


± 0.014 


pQCD 


+0.009 
-0.009 


+0.026 
-0.027 


+0.041 
-0.041 


+0.029 
-0.033 


Pow. corr. 


+0.060 
-0.040 


+0.055 
-0.037 


+0.081 
-0.054 


+0.125 
-0.084 


Exp. syst. 


+0.047 
-0.042 


+0.052 
-0.013 


+0.037 
-0.012 


+0.011 
-0.034 


Total 


+0.077 
-0.060 


+0.081 
-0.049 


+0.099 
-0.071 


+0.130 
-0.097 


ln(R) 


+0.003 


+0.019 


+0.023 


+0.012 


ln(R) mod. 


-0.005 


-0.015 


-0.025 


-0.019 


R 


+0.007 


+0.016 


+0.033 


+0.021 


R mod. 


-0.005 


-0.020 


-0.032 


-0.014 


x^ = 0.5 


-0.006 


-0.012 


-0.008 


-0.023 


XfI = 2.0 


+0.005 


+0.008 


+0.007 


+0.017 


M - 20% 


+0.060 


+0.055 


+0.081 


+0.125 


M + 20% 


-0.040 


-0.037 


-0.054 


-0.084 


91 GcV = OPAL 


+0.018 


+0.027 


+0.022 


+0.010 


91 GeV = ALEPH 


+0.018 


+0.020 


— 


— 


91 GeV = DELPHI 


< 0.001 


-0.013 


-0.010 


-0.017 


91 GeV = L3 


+0.021 


+0.024 


— 


— 


91 GeV = SLD 


+0.020 


+0.025 


+0.006 


-0.024 


y/s + 91 GeV 


+0.020 


+0.021 


+0.008 


-0.021 


v/i > 91 GeV 


-0.029 


-0.003 


-0.005 


+0.005 


y/s < 91 GcV 


+0.020 


+0.019 


+0.004 


+0.001 


fit range 


+0.018 
-0.031 


+0.008 
-0.001 


+0.030 
-0.005 


+0.002 
-0.017 



Table 2: Values of c*o (/•*/) derived using combined power corrections plus resummed QCD 
predictions to the event shape distributions of trust T, C-parameter, total and wide jet 
broadening Bt and Bw- The standard result for each variable is defined as the average 
of the results from the different matching schemes (ln(R), ln(R) mod., R, R mod.) using 
Xfj, — 1 and hi =2 GeV and considering the entire data set. In addition, the statistical and 
systematic uncertainties are given. Signed values indicate the direction in which a s (M z o) 
changed with respect to the standard analysis. Theoretical uncertainties are given by the 
matching scheme ambiguity, the choice of x^ and the uncertainty of the Milan factor M.. 
Experimental uncertainties are taken into account by varying the input data of the fit. All 
errors are treated as asymmetric uncertainties on a . 
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1-T 


C 


Bx 


B w 


ln(R) 
corr. 


283.7/277 
-16% 


162.9/170 
-60% 


161.0/171 
-11% 


162.1/132 

-51% 


ln(R) mod. 
corr. 


341.8/277 
-18% 


185.9/170 

-12% 


158.2/171 
-84% 


133.0/132 

-63% 


R 

corr. 


269.7/277 
-15% 


155.9/163 

-61% 


179.5/171 
-16% 


177.1/132 
-56% 


R mod. 
corr. 


329.0/277 
-75% 


171.4/163 
-15% 


160.7/171 

-84% 


139.3/132 

-64% 



Table 3: Values for x 2 /d.o.f. from the (a s ,o;o)-fits of combined power corrections plus 
resummed QCD predictions using the for different matching schemes. Also shown are the 
correlation coefficients of the fits. In the case of some C-parameter distributions, the fit 
ranges for the R- and the modified R-matching had to be chosen slighty different from 
those of the ln(R)- and the modified ln(R) matching in order to ensure the convergence of 
the fit. 
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Figure 2: Scaled distributions for the C-parameter measured by several experiments 
at ^/i = 35 to 183 GeV. The error bars indicate the statistical and experimental errors 
of the data points. The curves are the result of the simultaneous global (a s ,a;o)-fit using 
resummed QCD predictions with the modified ln(R)-matching plus power corrections. The 
fit ranges which are indicated by the solid lines are chosen individually for each center-of- 
mass energy. 
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Figure 3: Scaled distributions for Bt (left) and By/ (right) measured by several ex- 
periments at ^/s = 35 to 183 GeV. The error bars indicate the statistical experimental 
errors of the data points. The curves are the result of the simultaneous global (a s ,ao) - ht 
using resummed QCD predictions with the modified ln(R)-matching plus power correc- 
tions which include the revised power corrections to jet broadening observables [8]. The fit 
ranges which are indicated by the solid lines are chosen individually for each center-of-mass 
energy. 
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Figure 4: Results for a s (M z o) and a (2 GeV) from fits of power corrections to the event 
shape distributions of thrust T, C-parameter, total and wide jet broadening, B T and B w , 
respectively. The experimental and statistical uncertainties are represented by the solid 
error bars. The dashed error bars show the total error including theoretical uncertainties 
from the resummed QCD and the power corrections. The shaded band shows the one 
standard deviation around the weighted average. 
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